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Abstract 

We derive a general and closed-form result for the success probability in downlink multiple-antenna 
(MIMO) heterogeneous cellular networks (HetNets), utilizing a novel Toeplitz matrix representation. 
This main result, which is equivalently the signal-to-interference ratio (SIR) distribution, includes 
multiuser MIMO, single-user MIMO and per-tier biasing for K different tiers of randomly placed 
base stations (BSs), assuming zero-forcing precoding and perfect channel state information. The large 
SIR limit of this result admits a simple closed form that is accurate at moderate SIRs, e.g., above 5 
dB. These results reveal that the SIR-invariance property of SISO HetNets does not hold for MIMO 
HetNets; instead the success probability may decrease as the network density increases. We prove that 
the maximum success probability is achieved by activating only one tier of BSs, while the maximum 
area spectral efficiency (ASE) is achieved by activating all the BSs. This reveals a unique tradeoff 
between the ASE and link reliability in multiuser MIMO HetNets. To achieve the maximum ASE while 
guaranteeing a certain link reliability, we develop efficient algorithms to find the optimal BS densities. It 
is shown that as the link reliability requirement increases, more BSs and more tiers should be deactivated. 
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I. Introduction 

Network densification via small cells - resulting in heterogeneous cellular networks (HetNets) 

- along with novel multiple antenna technologies are two of the key methods for achieving a 
lOOOx capacity increase in 5G networks El- By densely deploying different types of small cells 
overlaid on the existing macrocell network, spatial reuse and thus area spectral efficiency (ASE) 
can be significantly improved and uniform coverage can be provided Q-fiSl. Independently, 
multiple antenna technology, known as MIMO, is playing an increasingly important role in 4G 
cellular networks, and is expected to be even more essential for meeting 5G target data rates. One 
key such technique is multiuser MIMO, which allows a base station (BS) with many antennas 
to communicate simultaneously with numerous mobile units each with a very small number of 
antennas. Although multiuser MIMO - also known as space division multiple access (SDMA) 

- has been known for quite some time and previous implementation efforts have been relatively 
disappointing, enthusiasm has been recently renewed, as seen in the extensive recent literature 
on “massive MIMO” [0, as well as the very recent 3GPP standardization of full-dimension (FD) 
MIMO, which can support 64 antennas in a 2D array at the BS to communicate simultaneously 
with 32 mobile terminals 0- 

A. Related Work 

A key result from early HetNet analysis was the derivation of the signal-to-interference-plus- 
noise (SINR) distribution, also known as the coverage or outage probability, where the HetNet 
is characterized by randomly placed base stations forming K tiers, each tier distinguished by 
a unique transmit power and density, i.e. the average number of BSs per unit area |f£fl, |[9l. A 
key resulting observation was that the signal-to-interference ratio (SIR) distribution is invariant 
to the BS densities, as long as the mobile connects to the BS providing the strongest received 
signal power. This SIR invariance property means that cell densification does not degrade the link 
reliability, and so the area spectral efficiency (ASE) of the network can be increased indefinitely 
by deploying more BSs. These early papers have resulted in a flurry of follow on work, e.g. 
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ma-ca, and see ffT6l for a survey. 

Multi-antenna transmissions bring significant additional complexity to HetNet analysis, pri¬ 
marily due to the complexity of the random matrix channel. As shown in iflTil . the invariance 
property may be lost in multi-antenna HetNets, i.e., the outage probability will increase as the 
BS density increases. This is mainly because the distributions of both the signal and interference 
depend on the number of BS antennas and the adopted multi-antenna transmission strategy 
of each BS. The work Q3 relied on stochastic ordering to compare different multi-antenna 
techniques, but such a method cannot be used for quantitative analysis, since the SINR and SIR 
distributions were not provided. That work was extended to incorporate load balancing and thus 
the achievable rate in lfl8ll . Other notable efforts on MIMO HetNets include work limited to two 
tiers [|T9l . [|20ll . and the analysis in ETTl . which focused on the interference distribution. In this 
paper, we attempt to solve a more fundamental problem, namely the determination of the SIR 
distribution in a MIMO HetNet, particularly a multiuser MIMO HetNet, and the resulting effect 
of the BS density on the network performance. By “performance”, we are interested in both link 
reliability - i.e. the success probability - and the sum throughput, which we characterize per 
unit area as the ASE. 

The link reliability vs. ASE tradeoff discussed in this paper is related to the notion of 
“transmission capacity” in wireless ad hoc networks li22l . If23l . with related multi-antenna results 
such as It24l - ll271 . In wireless ad hoc networks, to maximize the spatial throughput, i.e. ASE, 
while guaranteeing the link reliability at a certain SINR, the density of transmitters cannot 
exceed a certain value, which is called the transmission capacity lf23l . Naturally, allowing more 
simultaneous transmitters will increase the spatial reuse efficiency, but the interference to the 
receivers will become higher and so the SINR and thus link reliability will decrease, similar to 
the ASE vs. link reliability tradeoff studied in this paper. 

B. Contributions 


The main contributions of this paper are summarized as follows. 
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1) Tractable Expressions of the Success Probability for Downlink Multiuser MIMO HetNets: 
An exact expression of the success probability of the typical user in a general A'-ticr multiuser 
MIMO HetNet with zero-forcing (ZF) precoding and perfect channel state information (CSI) is 
derived. This expression has a simple structure and can be utilized for performance evaluation 
and system design of general MIMO HetNets. Moreover, based on this exact expression, an 
asymptotic (large SIR) expression of the success probability is also provided, which has a simple 
and symmetric form with respect to the network parameters. 

2) Key Properties of the Success Probability: Based on these tractable expressions, we ex¬ 
plicitly show the effect of the BS densities on the success probability, i.e., the link reliability. 
First, we prove that increasing the BS density of one tier will either increase or decrease the 
link reliability. Thus, the SIR invariance property no longer holds in MIMO HetNets. Second, 
we prove that the maximum link reliability of the network is achieved by only activating one 
tier of BSs. This is a surprising result, as nowadays a common intuition is that deploying more 
and more BSs will improve the network performance. Moreover, we show that such phenomena 
result from the fact that different tiers in MIMO HetNets have different influences on the link 
reliability, and the overall link reliability is an average of users in different tiers. 

3) ASE vs. Link Reliability Tradeoff: Since adjusting the BS densities will affect both the 
link reliability and ASE, a tradeoff is needed to balance these two aspects. Thus, we propose 
efficient algorithms to find the optimal BS densities that can achieve the maximum ASE while 
guaranteeing a certain link reliability requirement. If the required link reliability is low, all the 
BSs should be activated to achieve the maximum ASE. On the contrary, if the required link 
reliability becomes higher, more BSs should be deactivated. The extreme case is that to achieve 
the highest link reliability, only a single tier of BSs should be active. Therefore, care is needed 
when densifying a HetNet with multi-antenna BSs. 
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C. Paper Organization and Notation 

The rest of the paper is organized as follows. Section |TT] presents the system model. Section HIH 
derives the expressions of the success probability. In Section HV] we investigate the tradeoff 
between the ASE and link reliability. Finally, Section [V] concludes the paper. 

In this paper, vectors and matrices are denoted by hold-face lower-case and upper-case letters, 
respectively. The M x M identity matrix is denoted by 1 M . We use (•) 7 and (•)* to denote 
transpose and conjugate transpose, respectively. The Euclidean norm of the vector x is ||x||, 
while ||A||is the L\ induced matrix norm, i.e., ||A|| x = maxKj<„ Y^iLi l a *il f° r A £ R mxn . 
The expectation is denoted as E [•]. The notation X ~ Y means that X is distributed as Y. The 
normalized sine function is defined by sine (5) = . Finally, O (•) is the Landau’s notation, 

where f (x) = O (g (x)) as x —* x 0 means that there are constants C > 0 and £ > 0, such that 
for 0 < \x — Xq | < £, |/ (x)| < C \g (x)|. The asymptotic equivalence is denoted as ~, where 
/O) ~ g (x) as x —> oo means that f (x) / g (x) tends to unity as x —> oo. 

II. System Model 

In this section, we will first introduce a stochastic model for downlink multiuser MIMO 
HetNets, and then the main performance metrics considered in the paper will be presented. 

A. A Tractable Model for Downlink MIMO HetNets 

We consider a downlink cellular network consisting of K different tiers of BSs, indexed by 
K, — {1, 2,..., K}. In the /. th tier, the BSs are spatially distributed according to a homogeneous 
Poisson point process (PPP) |[28ll in M 2 of density X k , denoted by Each BS in the fcth tier 
has transmit power If, with M k antennas, and serves U k (U k < M k ) users at each time slot, 
i.e., intra-cell SDMA is considered. Note that when U k = 1, V/c £ /C, the network becomes a 
TDMA HetNet. Meanwhile, we assume that mobile users are distributed as a homogeneous PPP 
<f> u of density \ u , which is independent of 'Tfc, Mk £ /C. Each user has a single receive antenna, 
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Fig. 1. A demonstration of a 3-tier HetNet, in which each macro BS has 8 antennas and serves 4 users, each micro BS has 4 
antennas and serves 2 users, and each pico BS has a single antenna and serves 1 user. 


5S0. 


and is associated with one BSp]. We assume the network is fully-loaded (A u A k ,\/k G 1C), 
i.e., there are at least U k users in each cell in the /. th tier, and each BS will randomly choose 
Uk users to serve at each time slot. Note that this is a common assumption in the analysis of 
random cellular networks m, eel and the analysis of the non-fully-loaded network is left to 
our future work. An example of a 3-tier multiuser MIMO HetNet is shown in Fig. Q] 

The network is open access, which means a user is allowed to access any BS in any tier. 
Particularly, each user will listen to the pilot signals from different BSs, and measure the long¬ 
term average received power. We will consider a commonly adopted cell association rule [®, 
i.e., the user will be associated with the kth tier it; 


k = arg max RTFr,- ", 
j&K J J J 


( 1 ) 


where r 3 denotes the minimal distance from a user to its nearest BS in the jth tier, a > 2 is the 


'Note that our analysis also applies to single user MIMO, where each user in the fcth tier has Uk uncorrelated receive antennas, 
and BSs in the kth tier apply equal power allocation to the Uk streams. Then the success probability to the typical user analyzed 
in this paper becomes the success probability per stream. 

2 In the user association procedure, the first antenna normally uses the total transmission power of a BS to transmit reference 
signals for biased received power determination according to the LTE standard HD- 













6 


A 


path loss exponent and Bj is the bias factor of the j th tier, which is used for load balancing ||8]| , 
[|T6l . It is shown in J8| that the typical user will be associated with the kth tier with probability 

\PlB‘ 


Ai- = 


A 


( 2 ) 


where 8 = - and A = Y't-i A ;P?Bf. 

For each BS, equal power allocation is assumed in this paper, i.e., the BS in the /. th tier will 
equally allocate Pi- to its Uk users. Furthermore, we adopt ZF precoding and assume perfect CSI 
at each BSD Therefore, the received power of the typical user at origin o from a BS located at 
x G 'Ffc is given by 

Pi. 

, (3) 


P(x) = y^g x ,k IM 
Uk 


where g Xtk is the channel gain, and its distribution depends on whether the BS is the home BS 
(i.e., the associated BS) or the interfering BS. Assuming Rayleigh fading channels, it is shown 
in ll30l that g xk is gamma distributed with the shape parameter M k — U k + l if the BS at x is the 
home BS, i.e., g Xjk ~ Gamma (D k , 1), where D k = M k — U k + 1, while it can be approximated 
as g x , k ~ Gamma (U k , 1) if the BS is the interfering BS in the kth tier. This approximation has 
been commonly adopted in performance analysis of multiuser MIMO networks, e.g., |[T7l . llT9l . 
[f27l . m . In addition, the gamma distribution is a good approximation for the channel gain 
distributions for a large class of multi-antenna transmission schemes lf2TI . and thus the analysis 
in this paper can be extended to more general settings. 

In this paper, we assume that universal frequency reuse is applied, which implies that the user 
will not only receive the information signal from its home BS, but also suffer interference from 
all the other BSs. The resulting SIR of the typical user, served by the BS at xq in the kVa tier, 


’Note that in this paper, we assume the same path loss exponent in the network. However, the analysis can be extended to 
the model in which different tiers have different path loss exponents (e.g., ID). 

4 ZF is commonly used as the transmission strategy for multiuser MIMO in both the literature and industry, due to its simplicity 
and fairly well performance. It will be interesting to investigate how the findings in this paper will change if other precoding 
methods are used. However, this is left to our future work. 
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is given by 


SIR fc 



(4) 


Note that we consider SIR, rather than SINR, in this paper, as the SIR distribution captures the 
key tradeoffs in reasonably dense deployments. 

Remark 1: In summary, the most important parameters of the network are the per-tier pa¬ 
rameters: {Afc, Pfc, Bk, Mk, Uk}. Compared with existing works, the if-tier HetNet model in this 
paper includes intra-cell multi-antenna transmissions, which generalizes the SISO HetNet model 
(where M k — U k — 1 for all k e 1C) [[ 8 ft. Based on this model, the effects of the multiple antennas 
and the number of served users can be evaluated. Meanwhile, the derivation and analysis will 
be more difficult, as both the signal and interfering channel gains become gamma distributed, 
while in SISO HetNets, all the channel gains are exponentially distributed. 

B. Performance Metrics 

In this paper, we consider two performance metrics. The first one is the success probability, also 
called the coverage probability [[81, [[32ft . which measures the reliability of the typical transmission 
link. The definition of the success probability for the typical user is given by 


Ps = P (SIR > 7 ), 


(5) 


where 7 is the SIR threshold. Since the success probability is different if the typical user is 
associated with different tiers, the overall success probability in ([5]) can be written as in [H as 


K 



( 6 ) 


fc=i 


where p s (k) is the success probability when the user is served by the BS in the Ath tier. 


On the other hand, to describe the network capacity, we use the ASE as the mctriqj. defined 


by m, DU, El as 


K 



(7) 


k =1 


^Note that in some literature, e.g., ED and ED, the Shannon throughput of a cell, defined by E [log 2 (1 + SIR)], is used as 
the performance metric. The extension of our results to this case will be left to future work. 



for which the unit is bit/s/Hz/m 2 . Note that ASE measures the total data rate in a unit area 


normalized by the bandwidth, and it can directly tell the capacity gain by cell densification. 
The ASE defined in ([7]) applies to the fully-loaded networks. Note that if the network is not 
fully-loaded, the inactive BSs will increase the SIR and thus the success probability, but decrease 
the achieved ASE since there are fewer active links per unit area lfl3l . 

The success probability and ASE are two different aspects of a communication system, and 
both of them are fundamental metrics. A higher ASE means the network can support more users, 
i.e., with a higher spatial reuse efficiency, while a higher success probability implies that the 
transmission links are more reliable, i.e., a better quality of experience (QoE). While previous 
studies revealed the invariance property of the success probability in SISO HetNets, in this paper, 
we will investigate the interplay of these two metrics in more general multiuser MIMO HetNets. 
Since the ASE depends on p s (/;:), we will first derive it in Section [TTT1 and then investigate the 
tradeoff between ASE and p s in Section [TV] 


III. Analysis of the Success Probability 


In this section, we will first derive an exact expression of the success probability. Then, an 
asymptotic expression will be provided, which will expose its key properties in MIMO HetNets. 

A. The Exact Expression of p s 

To obtain the overall success probability p s from ©, we need to first derive the per-tier values 
p s ( k ). Assuming the typical user is served by the BS in the kth tier, then based on ©, p s ( k ) 
is given by 



( 8 ) 


The main difficulty of the derivation is the gamma distributed channel gains g X(j ± and g x j. There 
exist some previous works studying similar network models. While only a qualitative comparison 
between SDMA and TDMA was considered in ifTTl . the complicated expression of p s derived in 



[Tl8l based on the Faa di Bruno lemma obscured useful design insights. In this paper, we adopt 
a novel method and derive a more tractable expression for p s , as shown in Theorem Q] 

Theorem 1: The success probability of the typical user served by the BS in the fcth tier is 
given by 

Ps{k) = ^4 HQ^illi , ( 9 ) 


where A = Ylf=i $ = «’ D k — — 14 + 1, 11*11 x is the Li induced matrix norm, and 

Qo k is the D k x D k lower triangular Toeplitz matrix, given by 


in which 




Qo 

Qi % 

<?2 qi qo 

qn k -1 ••• qi qo 


( 10 ) 


qi = 


Va p s p s r ( Uj + ^ 5 

j-t J j j r (Uj)r(i + i)5-i 


U k B k A \ 

J^bP) 


r iB\ 


— 5, Uj + i; i + 1 — S; 


U k B k \ 

uiBp)' 

( 11 ) 


where 2 F± ( a,b;c;z ) is the Gauss hypergeometric function J35l P- 1005]. 

Proof: See Appendix [A] ■ 

Note that expression ([9]) is similar to the expression of p s of multi-antenna single-tier networks 
lf36l . However, compared with the result in lf36ll . the expression in Theorem Q] is in a simpler 
form while dealing with more general cases. Furthermore, the derivation provided in Appendix lA] 
adopts a different method which is more tractable than the one in [f36ll . Accordingly, the success 
probability of the typical user is given by 

K 

Ps — ^kP k F k 11 Q D \ 11 1 • ( 12 ) 

k =1 

The success probability of the single-tier network is given as a special case in the following 


corollary. 







10 


Corollary 1 : The success probability of the typical user in the single-tier network where each 
BS has M antennas and serves U users is given by 


Ps (k) = 


Qn 


( 13 ) 


where D — M — U + 1 , and Q d has the same structure as Qn k in (fTOl) . with the elements q, 
given by 


* = mwh ^ Fl +1 - * ■ 04 ) 

It shows in Corollary [T| that the success probability of the typical user is independent to the 
BS density, which indicates that the SIR-invariance property still holds in SDMA homogeneous 
cellular networks, which is consistent with SISO networks Il 32 l . 

Theorem Q] provides a tractable expression that can be easily evaluated numerically. However, 
it is still in a complicated form, which makes it difficult to directly observe the effects of 
different system parameters. To overcome this difficulty, in the next subsection, we will provide 
an asymptotic result of p s as the SIR threshold becomes large. 


B. The Asymptotic Expression of p s and Its Properties 


In this subsection, we will first derive an asymptotic expression of p s as the SIR threshold is 
large, i.e., 7 —* 00. Then, we will provide some basic properties of p s , which can help to better 
understand multiuser MIMO HetNets. 


The Taylor expansion of the coefficients q, in (fill) gives 

k /TT ^ x <5 


3 = 1 


U k B k \° 5 r(i + l- 5 )T(Cfj + 5 ) 


UjBj 


7 




( 15 ) 


s-i r (i + 1) r (Uj) \7 

Then, under the assumption that 7 —> 00, the asymptotic expression of the success probability 
is given in the following theorem. 

Theorem 2 : The asymptotic expression of p s (. k ) as 7 —> 00 is given by 

dDk+s'. 

r (D k ) 


Ps(k) rsj Ay d sinc ( 5 ) 


(U k B k )~ s 


EL aW g 


r(Uj+s ) 

mA 


( 16 ) 
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Fig. 2. The success probability with different SIR threshold, with a = 4, [Ai,A 2 ,A 3 ] = [1,5,10] x 10“ per km 2 , 
[Mi, Afe, Jtf 3 ] = [8, 4,1], [Ui,U 2 , U 3 ] = [4, 2,1], [Pi, P 2 , P 3 ] = [6.3, 0.13, 0.05]W and B k = 1 /U k . 


Proof: See Appendix [B] 

Based on Theorem [21 the success probability of the typical user is given by 

a ( p A 5 T ^ Dk+5) 

2^k=i A ky Uk J r (D k ) 


Ps = 


^ A k Ps (k) ~ 7 6 sinc (5) 


k =1 


Ef=l Ay 


r (jjj+s) 


(17) 


'3=i.-'3\Uj) T{Uj) 

In Fig. [21 we compare the asymptotic result in (fT71) with the exact numerical result in 
Theorem Q] as well as the simulation result. We consider a 3-tier HetNet, and the parameters 
are provided in the caption of the figure. From Fig. [21 we can find that the numerical result in 
Theorem Q] is exactly the same as the simulation result. Moreover, the asymptotic result in ([T71) 
provides as an upper bound when the SIR threshold 7 is not large, but the gap becomes smaller 
when 7 increases. This result confirms the effectiveness of the asymptotic result when the SIR 
threshold is large. For example, in this setting, the asymptotic result is almost the same as the 
exact result when 7 > 5 dB, and above 10 dB, is indistinguishable. Note that MIMO techniques 
will enhance the data rate at relatively high SIRs, so 7 > 5 dB is practical in MIMO HetNets. 
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From ( [T7T) . we can see a clear relation between the success probability and each parameter. It 
is a perfectly symmetric structure, since both the numerator and denominator are summations of 
each tier’s density {\k}_, the transmit power to each user |^|, and the effects of link array gain 
{D k } and multiplexing gain { U k }, respectively. Interestingly, this expression bears a similar form 
with the one in fU, which considered SISO HetNets. Based on (fT71) . we can easily investigate 
the effects of different parameters on the success probability in the general MIMO HetNets. 
Due to space limitations, in this paper, we will focus on how the BS densities {A/,} affect the 
success probability, and present two important properties. For convenience, we rewrite (fT71) in 
the vector form by defining the column vectors A = [Ai, A 2 ,..., A K ] T , c = [ci, c 2 ,..., c K ] r 
where q = 1 , and d = [d 1 ,d 2 ,... ,d K ] T where d t = ■ Then, (Q2]> is 

equivalent to 


p s ~ 7 ^sinc (5) 


c 1 A 
cFA' 


(18) 


1) Monotonicity of p s with respect to X k : Based on (fl8l) . we can obtain the following result. 

Lemma 1: When the SIR threshold 7 is large, if the ratio jj- = j s the same for 

all the tiers, p s is invariant with A/,,, V/c e /C; Otherwise, p s is monotonic with respect to A/,, i.e., 
p s will either increase or decrease as \ k increases. 

Proof: See Appendix O ■ 

This result explicitly shows that there is, in general, no invariance property for p s in MIMO 
HetNets. Deploying more BSs of one tier will either increase or decrease the success probability. 
The invariance property that p s is independent of the BS densities in SISO HetNets is only a 
special case, where & = 1 for all the tiers. Therefore, we should carefully consider how the 
success probability will be affected when densifying the HetNet with multi-antenna BSs. 

2) The maximum p s of the network: In this part, we will determine the maximum p s that can be 
achieved in a given MIMO HetNet. This is equivalent to considering the following optimization 


6 In this paper, we use {Afe} to denote the set {Xk : k G 1C}, representing the BS densities of all the tiers, while Xk is used 
to denote the BS density of the fcth tier. 
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problem: 


maximize p s (19) 

{•V} 

subject to 0 < A k < A™ ax ,Vfc <E /C, 


where A/, : can be regarded as the active BS density of the kth tier, and A“ ax is the actual BS 
density. The solution of this optimization problem is provided in Lemma [2} 


Lemma 2: The maximum p s with respect to {A^} in MIMO HetNets is given by 


pf ax = p s (fc) for k 


arg max 

3 


T (D, + 5) /r (£>,) 
r (u, + S) /r (u s ) ’ 


( 20 ) 


and the optimal BS density is 


I e (0. A?“] 
* 1 =° 

Proof: See Appendix [Dl 


U — arcr max ■ r (- D j+ < 5 )/ r ( D j) 

ft, — dig max, r(t/j+<5)/r([/,) ’ 


k arg maxj 


r(D 7 +«5)/r(D ; ) 

r ( u j+ s )/ r ( u j) ■ 


From Lemma [2l we see that when the numbers of antennas and the numbers of served users 


are determined, i.e., when { D k , If} are determined, then no matter how to change the BS 
densities, the success probability of the network cannot exceed the value in (l20l) . Furthermore, 

the maximum p s is achieved by only activating one tier of BSs, which has the largest value of 

r(A+£)/r(A) 

T(Ui+8)/Y(Ui) ' 

The two properties above have shown the effect of the BS densities on the success proba¬ 
bility in a general MIMO HetNet. In the next section, we will provide a more comprehensive 
investigation and jointly analyze the effects of the BS densities on the success probability and 
ASE. 


IV. ASE and Link Reliability Tradeoff 

In this section, we will investigate the effects of the BS densities on both the link reliability 
and ASE, and show there is a tradeoff between these two metrics. Moreover, efficient algorithms 
will be proposed to find the optimal BS densities. To start with, we will first consider a special 
case where all the BSs serve the same number of users. The purpose of investigating such a 
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special case is to reveal the physical insight of the tradeoff. Then, we will show the tradeoff in 
general MIMO HetNets using the asymptotic results provided in the last section. 

A. ASE and Link Reliability Tradeoff in U-SDMA Networks 

In this subsection, we will consider the tradeoff between the ASE and link reliability in the U- 
SDMA network lt37l . where each BS serves U (14 = U < M k for all k E /C) users. By limiting 
each BS to serve the same number of users, exact expressions of p s and ASE are available. 
Such a network still possesses the key characteristics of the MIMO HetNet. For example, when 
U — 1, the network becomes a TDMA HetNet. For simplicity, we assume B k — 1 for VA- e /C. 
Then p s (k) can be directly obtained from Theorem [H as given in the following corollary. 

Corollary 2: The success probability of the typical user served by the k\h tier in the unbiased 
(7-SDMA HetNet is given by 


Ps (k) — Q Dk , 


( 21 ) 


where Q Dk has the same structure as Q Dk in ([TO . with the elements q, given in (fT4l) . 
Therefore, the success probability of the typical user is given by 



( 22 ) 


and its vector form is given by 



(23) 


where the Ath elements of c and d are c k = p s (/;:) and d k = (jp) 5 , respectively. 
From Corollary O we have the following observations. 

Lemma 3: p s (/;:) in unbiased (7-SDMA HetNets has the following properties: 

• p s (k) is independent of the BS densities {A/,.} and the transmit power {P fc }. 

• Ps ( i ) > Ps ( j ) if and only if M, > Mj, and the equality holds only if M, = Mj. 


This result means that the tier with more antennas at each BS provides higher link reliability, 
since the interference suffered by each user has the same distribution in unbiased 1/-SDMA 
HetNets, while more transmit antennas will provide a higher diversity and array gain. While 
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p s (k) is independent of both {A/,} and { P k }, the overall success probability p s in (l22l) depends 
on {A k,Pk}, and p s can be regarded as a weighted sum of p s (k) for k € /C. From (1231) . the 
network success probability has the following properties: 

Lemma 4: p s in unbiased (7-SDMA HetNets has the following properties: 

• p s is mono tonic with respect to X k for V/c E 1C. 

• The maximum p s is achieved by only activating one tier of BSs which has the largest 
number of antennas, i.e., 

p“ ax = p s (k) for k = arg max Mj. (24) 

3 

Proof: Since (l23l) has the same structure with (fT8l) . Lemmas [H and [2] can be applied. 
Moreover, we have kf = p s (k) in this case, and from Lemma [3l we know the maximum p s ( k) 
is determined by M k . ■ 

By now, it is clear that the success probabilities in different tiers (p s (k)) are different, and 
p s is an average of each tier’s performance. Thus, the densification of different tiers will have 
different effects on p s . In particular, increasing the BS density of the tier with a lower p s (k) 
will pull down the overall p s . 

On the other hand, increasing the BS density will increase the ASE, which can be found in 
its expression: 

K 

ASE = U log 2 (1 + 7 ) ^ A kPs (k). (25) 

k= 1 

Note that since p s (k) is independent of the BS densities, (l25l) shows that activating as many 
BSs as possible can achieve a higher ASE. Combining the above results, we have the following 
two conflicting aspects: 

• To achieve the maximum ASE, it is optimal to activate all the BSs, but in this case p s may 
be low. 

• To achieve the maximum p s , activating only the tier with the largest number of antennas is 
optimal, but the ASE will be low. 
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Thus we need to investigate the tradeoff between ASE and link reliability. We are interested in 
maximizing the achievable ASE given a requirement on p s , formulated as the following problem. 


T 0 : maximize ASE 

(A 1 

subject to p s > 0, 

0 < A fc < Aj^.Vfc e /C, 


(26) 


where ASE is given in (l25l) and p s is given in (l22l) . 

Note that in unbiased [7-SDMA networks, Problem V„ is a linear programming problem [(3811 . 
which can be solved efficiently. To get more insights, we will derive a more explicit solution. 
From Lemma lU we know that if 0 > max/, p s (k), there is no feasible solution, which implies 
that whatever the BS density is, the network cannot achieve such a link reliability requirement. 
Thus, in the following analysis, we assume 0 < max fe p s (/c). 

Denote y k = p s (k) Afc, then Problem Vo is equivalent to 

Vu -sdma : maximize Ylk=i Vk (27) 

{yk} 

subject to Y.k=i b kVk > 0, 


0 <s/fc<j/r x ,Vfee/C, 


Pk\ s 


where b k = (^) 


1 - 


e 

Ps(fc) 

K 


. Note that b k < 0 if p s ( k ) < 0. To find the solution, we start 


from {y* = yf**}. If E i= i b jVj > then the optimal solution is to activate all the BSs. On 


the contrary, if Ylf= i b jVj < it means we need to decrease some y* k until b jVj = 0- 

Without loss of generality, assume some of {b k }, i.e., b n } where n < K, are negative, 

and b\ < b 2 < b n < 0. Then, we need to firstly decrease y\, as b\ has the most negative effect. If 


K 


3/i = 0 and b jV*j < 0^ then y\ should be decreased. So on and so forth until i b jUj = 0 
The formal solution is described in Algorithm [H 


K 


From Algorithm Q] we see that the value of b k is related to the transmit power P k and the 
number of antennas M k . A negative b k will have a negative effect on p s . If the BSs in the A th 
tier have a smaller M k and a larger P k , the value of b k will be negative and smaller. Therefore, 
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Algorithm 1 Finding the optimal solution of Problem V 0 in unbiased U -SDMA HetNets 
l: Initialize y\ <— y[ nax for k E 1C and n <— 1; 

2: while Ef=r hy* k < 0 do 

3: i <— the index of the nth minimal value among {b k }; 

4 : Vi <- 0 ; 

5: if Ef=i hyt > o then 

6: v* <- i Ef=i hy%\ 

7: break; 

8 : end if 

9: n <— n + 1; 

to: end while 


to achieve a high link reliability requirement, the BSs of the tiers with a small b k should be 
switched off. Note that one special case is that all the BSs have the same number of antennas, 
i.e., M k = M for \/k E 1C. In this special case, p s (k) is the same for all k E 1C. Thus, all the 
values of {b k } are either greater than 0 or less than 0, and there is no tradeoff since p s becomes 
a constant. SISO HetNets belong to this special case. 

B. A Demonstration of the Tradeoff in Unbiased U-SDMA Networks 

In this subsection, we will demonstrate how to achieve the maximum ASE given the link 
reliability requirement in a 3-tier HetNet, consisting of micro-BSs, pico-BSs and femto-BSs, 
where [Mi,M 2 ,M 3 ] = [4,2,2], [Pi , P 2 , Vo] = [6.3,0.13,0.05] Watts ||39l , and the actual BS 
densities are [A“ iax , A™ ax , A™ ax ] = [1,2,5] x 10 2 per km 2 . We consider two cases with U = 1 
and U = 2, respectively. 

In Fig. [3l we show the tradeoff between the maximum ASE and the link reliability requirement 
0. As a benchmark, we also consider the SISO HetNet with the same transmit power {P k } and 
actual BS densities {A™ ax }. From Fig. [3j we can find that: 1) Compared with the SISO HetNet, 
deploying multi-antenna BSs can increase both the ASE and link reliability. For example, in the 
SISO HetNet, p s is always 0.56 whatever the BS density is. But when using multi-antenna BSs, 
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Fig. 3. The maximum ASE with different requirements of the success probability, with a = 4, 7 = 0 dB. [Mi, M 2 , M 3 ] = 
[4,2,2], and [Pi,P 2 ,P 3 \ = [6.3,0.13,0.05] Watts. The actual BS densities are [A™ 311 , A™ ax , A™ ax ] = [1,2,5] x 10 2 per km 2 . 


the network can achieve better ASE and p s . 2) Comparing the cases with U — 1 and U = 2, we 
see that serving one user at each time slot (i.e., TDMA) can obtain a higher maximum p s . This 
is because the channel gains of the interference links are small when U is small. However, the 
maximum ASE is higher when U = 2, compared with U = 1, since there are more links in a 
unit area. 

The corresponding optimal BS densities of the case with U = 2 are provided in Fig. |U 
When 0 > 0.58, using {A™ ax } cannot satisfy the link reliability constraint. From calculation, we 


p k \ s 


find that the 2nd tier has the minimal value of (^r) 


1 - 


e 

Ps(fc) 


, which means we need to first 


decrease the BS density of the 2nd tier, as confirmed in Fig. [4] Comparing the 3rd tier with the 
2nd tier, we can find that M 2 = M 3 , but P 2 > Pi, which means that BSs in the 2nd tier have 
more negative effects on the link reliability as they will cause higher interference. This is the 
reason why we need to decrease the BS density of the 2nd tier to improve the link reliability. 
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Fig. 4. The optimal BS densities with different requirements of the success probability, with a = 4, 7 = 0 dB, [Mi, M 2 , M 3 ] = 
[4,2,2], U = 2, and [P\, P 2 , Pz\ = [6.3,0.13,0.05] Watts. The actual BS densities are [A™ 331 , A™ 8,31 , A™ ax ] = [1,2,5] x 10 2 
per km 2 . 


C. ASE and Link Reliability Tradeoff in General MIMO HetNets 


Now, we will concentrate on the tradeoff in a general multiuser MIMO HetNet. Since the 
analysis becomes extremely difficult with the exact expression of p s in ([9]). we resort to the 
asymptotic expression in (fT6l) . By substituting (fT6l) into Q), the ASE is given by 

(cjx) (4 A) 


ASE ~ 7 ^sinc (5) log 2 (1 + 7 ) 


d T A 


(28) 


where cq = [c n , C12, • • • Cik] T , and its element is given as cu = P^Bf, while c 2 = [c 2 i, c 2 2, ... c 2 k] T , 
in which c 2t = U }~ 5 and d has the same expression as in (fl8l) . 

From (l28l) . we can find that the ASE is a quadratic-over-linear function with respect to A. It 
means, in general, we cannot guarantee that increasing the BS density will always increase the 
ASE, and the effect of the BS density on the ASE becomes complicated. In the following, we 
consider the same problem V Q as in (l26l) . to find the optimal BS densities that can achieve the 














20 


maximum ASE. Based on (fT8l) and (l28l) . the optimization problem is given as 

V\ : maximize — — (29) 

subject to 7 _,5 sinc (5) > 0, 

0 < X k < Xf**, k e K. 

Since the maximum p s is obtained by activating only one tier which has the maximal value of 
r^+^/rltr) ( C L Lemma [2]), and p™ ax ~ 7 _<5 sinc (5) max, rffi+sj/rxrr) » we only consider the 
case when 0 < p™ ax . 

Problem V\ is a non-convex problem since the ASE is not a concave function with respect 
to the BS density. Fortunately, Dinkelbach has proposed an algorithm to solve the nonlinear 
fractional problems li40l . By defining N (A) = (c{ A) (c^A) and 

F (t) = max {N (A) - td T A | A e S} , (30) 

where S — | A | 7 _<5 sinc (5) > 0, 0 < A fc < A™ ax , k G /C j, the optimal BS density A* can be 

obtained by finding t* such that F (£*) = 0 flOl . Specifically, by iterating fW — _/y (A*) / (d 7 A*), 
where A* is the optimal solution of the right hand side of (l30l) . F will converge to 0 and 
A* will be the optimal BS density. However, the optimization problem in (l30l) is still a non- 
convex problem due to the non-convex function N (A). To resolve this difficulty, we resort to 
the sequential convex programming (SCP) by approximating N (A) with the first order Taylor 
expansion, given by N (A) « N ^A^j + VN ^A (n )j ^A — A^j. It is a simple but effective 
method and has wide applications, e.g., see tl4Tll . Therefore, given the nth iterative \ {v) and the 
zth iterative l (l K the convex optimization problem is given by 

V 2 (X (n \t (l) ^ : maximize N (a w ) + ViV (a ( n) ) T (a - A (n) ) - t«d T A 
subject to ( 7 _,5 sinc ( 5 ) c T — 0d T ) A > 0, 

0 < A < A max . 


Then, the optimal BS density can be obtained by Algorithm [2l Note that since SCP can only 
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obtain a local maximum, we will randomly generate multiple initial values of A' 0 - 1 to find a 
better solution. 


Algorithm 2 The Locally Optimal BS Densities in General MIMO HetNets 
l: Initialize A^ random value between [0, A max ], n <— 0, i 0 and assign e a small value; 

(cfAW)(c^A(”)) i 


2: tW 




d T A*") 


4 : if 


A (n) 


/ 


A (n) 


: Solve Problem P 2 (Aand obtain the optimal value A*; 


> e then 


5: n <— n + l, A^ <— A*, and go to Step El 

6 : else 

7: if (cf A*) (4A*) - t«d T A* > £ then 

8: i <— i + 1, and go to Step El 

9 : else 


10: return A*; 

it: end if 

12: end if 


D. A Demonstration of the Tradeoff in General MIMO HetNets 

In this subsection, we will use Algorithm El to evaluate the tradeoff between the ASE and 
the link reliability. We consider a 3-tier HetNet, where [M 1 ,M 2 ,M 3 ] = [8,4,1], [Ui, C/ 2 , Uf\ = 
[4,1,1], B k = 1/U k for k = 1,2,3, [Pi,P 2 ,P 3 ] = [6.3,0.13,0.05] Watts, and the actual BS 
densities are A max = [A™ ax , A™ ax , A™ 3 ^] = [1, 5,10] x 10 2 per km 2 . For AlgorithmEl we generate 
20 randomly initial values of A (0) , and set e = 10~ 6 . In Fig. [5J we show the tradeoff between 
the maximum ASE and the link reliability requirement 0, while the corresponding optimal BS 
densities are shown in Fig. [6] 

First, we find from both figures that the proposed algorithm can achieve almost the same 
performance as using the exhaustive search, while the proposed algorithm runs much faster. 
Second, from Fig. [5J we find that in general cases, there exists a tradeoff between the ASE and 
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Fig. 5. The maximum ASE with different requirements of the success probability, with a = 4, 7 = 5 dB, [Mi, M2, M3] = 
[8, 4,1], [U\, U2, U 3 ] = [4,1,1], Bk = 1/Uk for k = 1,2, 3, and [Pi, P2, P3] = [6.3, 0.13, 0.05] Watts. The actual BS densities 
are A max = [1, 5,10] x 10 2 per km 2 . 



Fig. 6 . The optimal BS densities with different requirements of the success probability, with a = 4, 7 = 5 dB, [Mi, M2, M3] = 
[ 8 , 4,1], [Pi, U2, U3] = [4,1,1], Bk = 1 /Uk for k = 1,2,3, and [Pi, P2, P3] = [6.3, 0.13, 0.05] Watts. The actual BS densities 
are A max = [1, 5,10] x 10 2 per km 2 . 
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link reliability, and the higher the link reliability requires, the lower the ASE can be achieved. 
More interestingly, from Fig. [6} we find that: 1) Even in general HetNets, the maximal ASE 
is achieved by activating all the BSs. 2) With the increasing requirement of the link reliability, 
the BS density will decrease from one tier to another tier. Moreover, only when the BS density 
decrease to 0, the BS density from another tier will start to decrease, which is the same as U- 
SDMA HetNets. Thus, we can infer that for a given 0, different tiers have different influences 
on the network, and there is an ordering of such influences. Recall that in [/-SDMA networks, 


the ordering can be explicitly obtained from {b k }, where b k = (jf)' 


1 - 


PsO) 


, and the tier 


with the minimal negative b k among k G K, has the most negative effect on the link reliability. 
However, in the general MIMO HetNets, so far we are unable to derive an exact expression to 
find such an ordering, since both signal and interference distributions are complicated. But with 
our algorithm, we can numerically find the effects of different tiers on the link reliability. For 
example, from Fig. [6] we can see that tier 1 (circle points) has the most negative effect on the 
link reliability. When 0 > 0.5, the BSs from tier 1 need to be deactivated to guarantee the link 
reliability. It is because the interference caused by tier 1 is large due to the high multiplexing 
gain U\ — 4. 


V. Conclusions 

In this paper, we developed a new set of analytical results for performance analysis of downlink 
multiuser MIMO HetNets, where multi-antenna BSs use SDMA (ZF precoding) to serve multiple 
users. Both exact and asymptotic expressions of the success probability were provided. We 
focused on the effect of the BS densities, and proved that there is no invariance property for the 
success probability in general MIMO HetNets. A unique tradeoff between the link reliability and 
the ASE was revealed. By using the proposed algorithms, we found the optimal BS densities 
to achieve the maximum ASE while guaranteeing a given link reliability requirement. It was 
shown that the maximum ASE of the network is achieved by activating all the BSs, while the 
maximum link reliability is achieved by activating only one tier of BSs. 
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The link reliability vs. ASE tradeoff analyzed in this paper provides a new perspective in 
designing multiuser MIMO HetNets, which is fundamentally different from SISO HetNets. 
Both the link reliability and ASE should be considered when evaluating different transmission 
techniques for HetNets. A basic network model is considered in this paper, while more general 
models require further investigation, such as multi-slope path loss models ll42l . together with 
more sophisticated techniques such as interference coordination |[43l and load balancing lTT8l . 


Appendix 

A. Proof of Theorem [7] 

The success probability in ([8]) can be rewritten as 


K 


PiUk 


Ps (fc) = p I g X0 ,k >r%^2 j2fl 9x ’ j 


(31) 


j=i xe\iq\{xo} 

Since g XQ>k ~ Gamma (D k , 1), the above equality can be expressed as 

~D k — 1 


p s (k) = E s 


E E ' 

n=0 


«L e -u 


n\ 


(32) 


where s = and I = £* i E aG ^\{*o} pm 9 ** Ik 


. Note that for a fixed s, E/ [e _si ] is 
the Laplace transform of /, denoted as Cj (s). Following the property of the Laplace transform, 
we have E j [. I n eT sI ] = (— l) n (s), where ( s ) i s the nth derivative of Cj (s). Then, we 

have 


p s (k) = E s 


D k —1 

E 

n=0 


(n) 


n\ 


cr («) 


(33) 


To derive p s (k) based on (l33l) . we start from the Laplace transform of / for a fixed s, given 


by 


K 


PiU k 


I 1 j u k I, || —q 

exp | —s 2_^ 2^ -pjr9xjM 

j= 1 x6$j\{x 0 } ^ J 


jCj (s) = E 


(34) 
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Since {'I' J : j e K,} and {g x j : x £ 'I > b ,j € JC} are independent, the above equality can be written 
as 


K 


Cl (s) — JjEg 


(a) 


3 =1 

n k ' 

j=i 


n e 

xe'i'jMxo} 

n 

xG'I'jVtxo} 


exp 


PjU k „ . 
-s— — q x j x 
PkUj 13 " 1 


P Ui ^ ~ Uj 

l + s Jw, M 


(35) 


where (a) follows from g x j ~ Gamma (Uj, 1) for x £ '\> r Then, the Laplace transform of I can 
be derived using the probability generating functional (PGFL) Il44ll . which is given as 


K 


Ci ( s ) = exp < —TV A j 

3 = 1 


P,u u ^ 


i I i i — 7 ^ /c — — 

1 _ 1 + s ~pJT j v 1 


dv } , 


(36) 


where ry is the minimal distance between the typical user and the BSs of the jth type. Based 
on (l36l) . the nth derivative of £/ (s) with respect to s can be written as the following recursive 
form 


4 n) w = 


Ti — 1 

£ 

2=0 


77, — 1 

i 


K 

3 =1 


(-i 


PjU h -£L 


PkUj 


V 2 


dv 


, , PjU k _2A J 

1 + s 7 tm v 2 


Uj+n—i 


cf (s) 


(37) 


. The 


where (Uj) _. is the Pochhammer symbol. 

Denote x n = ^ (—s) n (s) for n > 0, then it shows that p s ( k) = E. s J2n=o 1 x n 

expression of x 0 can be obtained directly from (1361) . while substituting (1371) to x n , we can get 
for n > 1, 


x n = 7r 


n —1 

£ 

i=0 


77, — X 


Xi 


n 


K 

3 =1 




PjU , fe „ “ 




V 2 


dv 


r 2 [n — i) 


PjU k —— 
1 + 2 




(38) 


In the expression of x n , the integration limits are from r| to infinity. As the typical user is 
associated with the BS in the fcth tier, it implies that P k B k r^ a > PjBjrJ a for j £ /C, which 
is equivalent to rj > (77^j r^, where A = £ Moreover, substituting s = A /r k to xo, the 
expression of x 0 is given by 


K 

x 0 = exp ^ —7r Aj 

7=1 


2 


, , . a PjU k _^~ Uj 
1 + r Y r kTrrr v 2 

* k Uj 


dv > , 


(39) 
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and for n > 1, 


n —1 

X n = 71 X 
i=0 


n — i 

i 

n 


K 

£ A ; 

i=i 




a PiUk —— 
2 


2 77 — 7 


-I . ^ syPn XJlr - — 

1 + r r r ktm v 2 


k^j 


Uj+n—i 


dv 


(40) 


Denote q t as (flTT) . then with some manipulation, the expression of x n can be simplified as 


x 0 = exp 


7r rt 


HBl 


(Qo ~ A) 


and for n > 1, 


nrt 




n— 1 

£ 


77 — 7 


PuBi *■—' n 

* * i =0 




(41) 


(42) 


To solve the recurrent relation of x n , we define two power series F (z) and G (z) as 

OO OO 

f(z) = Y 1 ^ n , G w = E XnZl1 - ^ 43) 

71=0 71=0 

Using the properties of F' (z) = J2n=o Mn^ 1 and F (z) G (z) = J2n=o (ELo Qn-iXi) z n , from 


, we can get the equality as 

zG' (z) = 


nr i 


Hsi 


zF' (z) G (z) 


Then, by solving the above differential equation, we have G (z) = c exp 
c is a constant to be determined. Since G (0) = x 0 and G (0) = cexp 
we get c = exp (^-A^A). Thus, G (z) is given by 


7T rt 


6 nt> - 


P k B k 


(44) 


F (z)), where 


pM 


qo), from (flTT) . 


G (z) = exp 


nri 




(F(z)-A) 


(45) 


Dk~ 1 


Recalling that p s (k) = E rfc Yz n = 


0 


, we define 


T ( z ) — ^r k [G (z)] — £ t n z n . 


(46) 


n=0 

From [[8l Lemma 3], the probability density function of r k is given by 

-^ 2 £?= iA 


fr k ( r ) = 2vrr 


A / p U \ 5 

E A ' 33 


Then we can get 


T (z) = E, 


Zk 


.7=1 


exp 


PkB k 


K , ( P J B j' 5 


-3= 1 J V P k B, 


nrl 


PlFtl 


(F (z) - A) 


A 


F z) 


(47) 


(48) 
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Therefore, the success probability is given by 

Dk — 1 Dk~ 1 

Pi 


D k -1 


(^) — ~ ^ ^ l 2=0_ 


n=0 


n=0 


n: 




1 d" 


—■' n! ch n I -P 1 (- 2 ) 

n =0 v v ' 


(49) 


2 = 0 


Finally, from 1(431 p. 14], it can be shown that the first D k coefficients of the power series 
is the first column of the matrix Q D *, where Qo k is given in (fTOl) . Thus, (|49|) is equivalent 
to ©, which completes the proof. 


B. Proof of Theorem |2] 

By substituting (fT51) into the power series F (z), we can obtain 


oo K 

F(z) = Ei‘ z, ~'£ x i p i B i 

i =0 j=l 

and it can be expressed as 

F(A~Y.hpB‘ J 

3=1 

Thus, the power series T (z) is given by 

A 


U k B k \ 

6 r (Uj + <5) 

VjBj ~‘) 

T(Uj) 

U k B k A 

6 r (Uj + 5) 

If, BP) 

r m 


5 T (i + 1 — 5) ^ 

A — i Ft i -I- 1 1 


»=0 


T (z) = 


■^r(i-s) (-*- z ) 


-5 


F ( z ) 


sp K \ pS tj 5 ( UkBk xA r(Fj+<5) 

Z^j=i j ^ Uj b.j i) r (Uj) 

Based on the above expression, the coefficient t n is given by 


tn = ~J {n) (z) |,=C 


^n! (^)n r(l-S) 


n\ 


sp K \ pS p<5 f U k B k ~\ T(Uj+S) 

^3=1 U B 3\ifb-^) -m-y 


As p s (k) = En=o 1 *n» and the summation Y‘n=o ‘ ft (<*)„ = 


D k — 1 1 /'ST'i _ r(g fc +<5) 


(50) 


(51) 


(52) 


(53) 


is given by 


Ps(k) 


A 


n= 0 n\ \ In r(l+<5)r(D fe ) 


r (Dk+6) 


, the success probability 


r(i-<5)r(i+5)r(D fe ) 


\ pS p5 (UkBk A r(Fj+i5) 
2-/ \ A J l j n .i \U:I1, ) r {Uj) 


( 54 ) 


Using the equality T (1 + S) T (1 — 5) = 


nS 


sin(7T(5) sinc(<5) 


, we can obtain (fT6l) . 
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C. Proof of Lemma [Tj 

We consider the function y = where A = [Ai, A 2 ,. .. \r] T , c = [ci, c 2 ,..., c K ] and 
d = [di, d 2 , ■ ■ ■, d K ). The partial derivative of / (A) with respect to A* is then given by 


dy_ 

dK 


Si =1 didjXj ( 2 J. 


di 


(d T \y 


(55) 


It shows that changing A, will not change the sign of J^-, i.e., y is monotonic with respect to 


Aj. Moreover, y is independent of A if ^ ^ for all i,j e {1,2,..., K}. 


D. Proof of Lemma \ 2] 

In this proof, we consider a more general case where the optimization problem is given by 

'Pw/cons : maximize ^ (56) 

subject to AA < b, 

A > 0, 

where AA < b represents arbitrary constraints, A is a n x K matrix, and b is a n x 1 vector 
with positive elements. We will find the optimal solution A* in the following derivation. 

First we consider the optimization problem without A A < b constraint, i.e., the problem is 
given by 

V w /o cons : maximize {pA (57) 

subject to A > 0. 

Without loss of generality, we assume ^ ^ > • • • > g-, then from (1551) . we can find that to 

maximize p s , \ K should be 0 since < 0. 

Then, repeating the same procedure, we can find \k-i = 0, \r -2 = 0, and A 2 = 0 
successively. Finally, the objective function p s is equal to with any Ai > 0, which is the 
solution of V w/ocons . 



29 


Secondly, we consider the optimization Problem V w / cons , and we want to prove that the 
optimal A* of Problem V w / OCO ns is also the optimal solution of problem Vw/cons- To do so, we 
need to prove 1) A* is feasible for V w / cons , and 2) the optimal solution of Vw/cons is A*. 

To prove the feasibility, we assume A*, and A| are optimal solutions of V w / OC ons , but A A* > b. 
Since b > 0, and Ag = k\\ for any k > 0. Then 3k > 0, so that AAg = fcAA* < b, i.e., there 
exists an optimal solution of V w / ocon s , which is feasible to V w / CO n S - 

Finally, since V w / con s has one more constraint than V w / OC ons , the solution of V w / CO ns should 
be the subset of the solution of V w / OCO ns and the maximum value of V w / con s will no greater 
than the maximum value of V w / OC ons- As we have proved that the optimal solution of V w / OC ons 
is feasible in V w / cons , we can obtain that the maximum value of V w / CO ns is equal to 

Now we come back to the general MIMO HetNet case. The original problem (fl9l) indicates 
A = I K and '-j- = Vtu+ 3 )/r(u) • Therefore, the maximum p s is obtained by only deploying one 
tier which has the maximal value of V- 

di 
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